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The MiniBooNE experiment at Fermilab reports results from an analysis of £" e appearance data 
from 11.27 x 10 20 protons on target in antineutrino mode, an increase of approximately a factor of 
two over the previously reported results. An event excess of 78.4 ± 28.5 events (2.8a) is observed 
in the energy range 200 < E c j E < 1250 MeV. If interpreted in a two-neutrino oscillation model, 
— > v e , the best oscillation fit to the excess has a probability of 66% while the background-only 
fit has a x 2 -P r °bability of 0.5% relative to the best fit. The data are consistent with antineutrino 
oscillations in the 0.01 < Am 2 < 1.0 eV 2 range and have some overlap with the evidence for 
antineutrino oscillations from the Liquid Scintillator Neutrino Detector (LSND). All of the major 
backgrounds are constrained by in-situ event measurements so non-oscillation explanations would 
need to invoke new anomalous background processes. The neutrino mode running also shows an 
excess at low energy of 162.0 ± 47.8 events (3.4<r) but the energy distribution of the excess is 
marginally compatible with a simple two neutrino oscillation formalism. Expanded models with 
several sterile neutrinos can reduce the incompatibility by allowing for CP violating effects between 
neutrino and antineutrino oscillations. 



There is growing evidence for short-baseline neutrino 
anomalies occurring at an L/ E v ~ 1 m/MeV, where E v 
is the neutrino energy and L is the distance that the neu- 
trino travelled before detection. These anomalies include 
the excess of events observed by the LSND [T] and Mini- 
BooNE [2H1] experiments and the deficit of events ob- 
served by reactor [3] and radioactive-source experiments 
[6j. There have been several attempts to interpret these 
anomalies in terms of 3+N neutrino oscillation models 
involving three active neutrinos and N additional sterile 
neutrinos [7HT2]. (Other explanations include, for ex- 
ample, Lorentz violation |13) and sterile neutrino decay 
[14].) A main goal of MiniBooNE was to confirm or re- 
fute the evidence for neutrino oscillations from LSND. 
Of particular importance is the MiniBooNE search for 



— > v e oscillations since this was the channel where 
LSND observed an apparent signal. This paper presents 
improved results and an oscillation analysis of the Mini- 
BooNE v e appearance data, corresponding to 11.27x 10 20 
POT in antineutrino mode, which is approximately twice 
the antineutrino data reported previously [1]. 

Even though the first goal of this article is a presen- 
tation of the improved antineutrino results, a secondary 
goal is to contrast and compare these results with im- 
proved MiniBooNE neutrino measurements and, there- 
fore, the details of both the neutrino and antineutrino 
analysis will be given. Since the original neutrino result 
publication [3], improvements to the analysis have been 
made that affect both the v e and v e appearance search. 
These improvements are described and used in the anal- 



2 



yses presented here. 

The neutrino (antineutrino) flux is produced by 8 GeV 
protons from the Fermilab Booster interacting on a beryl- 
lium target inside a magnetic focusing horn set at posi- 
tive (negative) polarity. In neutrino (antineutrino) mode, 
positively (negatively) charged mesons produced in p-Be 
interactions are focused in the forward direction and sub- 
sequently decay primarily into (y^). The flux of neu- 
trinos and antineutrinos of all flavors is simulated using 
information from external measurements [15] . In neu- 
trino mode, the v^, v^, v e , and v e flux contributions at 
the detector are 93.5%, 5.9%, 0.5%, and 0.1%, respec- 
tively. In antineutrino mode, the v^, v^, v ei and v e flux 
contributions at the detector are 83.7%, 15.7%, 0.4%, 
and 0.2%, respectively. The and fluxes peak at 
approximately 600 MeV and 400 MeV, respectively. 

The MiniBooNE detector is described in detail in ref- 
erence [TBJ- The detector is located 541 m from the 
beryllium target and consists of a 40-foot diameter sphere 
filled with 806 tons of pure mineral oil (CH2). Neutrino 
interactions in the detector produce charged particles 
(electrons, muons, protons, pions, and kaons) which in 
turn produce scintillation and Cherenkov light detected 
by the 1520 8-inch photomultiplier tubes (PMTs) that 
line the interior of the detector and an optically iso- 
lated outer veto region. Event reconstruction and par- 
ticle identification are derived from the hit PMT charge 
and time information. In particular, the reconstructed 
neutrino energy, E® E , uses the measured energy and an- 
gle of the outgoing muon or electron assuming charged- 
current quasi-elastic (CCQE) kinematics for the event. 

The signature of —¥ v e and — > v e oscillations is 
an excess of v e and i/ e -induced CCQE events. Recon- 
struction |17j and selection requirements of these events 
are almost identical to those from previous analyses 0] 
with an average reconstruction efficiency of ~ 10 — 15% 
for events generated over the entire volume of the detec- 
tor. Recent improvements to the analysis include a bet- 
ter determination of the intrinsic v e background from K + 
decay through the measurement of high-energy neutrino 
events in the SciBooNE experiment [18], a better deter- 
mination of NC 7T° and external event backgrounds in an- 
tineutrino mode due to the increase in statistics of the an- 
tineutrino mode data sample, and the use of a likelihood 
fit with frequentist corrections from fake data studies for 
both the neutrino-mode and antineutrino-mode analy- 
ses. The detector cannot distinguish between neutrino 
and antineutrino interactions on an event-by-event basis. 
However, the fraction of CCQE events in antineutrino 
(neutrino) mode that are due to wrong-sign neutrino (an- 
tineutrino) events was determined from the angular dis- 
tributions of muons created in CCQE interactions and 
by measuring CC single tt + events [19]. 

The predicted v e and v e CCQE background events 
for the neutrino oscillation energy range 200 < E® E < 
1250 MeV are shown in Table U for both neutrino mode 



TABLE I: The expected (unconstrained) number of events for 
the 200 < E® E < 1250 MeV neutrino oscillation energy range 
from all of the backgrounds in the v e and v e appearance anal- 
ysis and for an example 0.26% oscillation probability averaged 
over neutrino energy for both neutrino mode and antineutrino 
mode. The table also shows the diagonal-element systematic 
uncertainties whose effects become substantially reduced in the 
oscillation fits when correlations between energy bins and be- 
tween the electron and muon neutrino events are included. 





Process 


Neutrino Mode Antineutrino Mode 




v„ & CCQE 


37.1 ± 9.7 


12.9 ± 4.3 




NC 7T° 


252.3 ± 32.9 


112.3 ± 11.5 




NC A — > N-y 


86.8 ±12.1 


34.7 ± 5.4 




External Events 


35.3 ± 5.5 


15.3 ± 2.8 




Other i>ft & 


45.1 ± 11.5 


22.3 ± 3.5 


V e 


& Ve from fj, Decay 


214.0 ± 50.4 


91.4 ± 27.6 


Ve 


& Ve from K Decay 


96.7 ± 21.1 


51.2 ± 11.0 


v e 


& v e from Kl Decay 


27.4 ± 10.3 


51.4 ± 18.0 




Other v e & v e 


3.0 ± 1.6 


6.7 ± 6.0 




Total Background 


797.7 


398.2 




0.26% V„ -> Ve 


233.0 


100.0 



and antineutrino mode. MiniBooNE does not have the 
electron versus gamma particle identification capabili- 
ties to determine whether observed events are due to 
charged-current (CC) electron events, as expected for an 
oscillation signal or intrinsic beam v e /v e background, or 
to background gamma events from neutral-current (NC) 
interactions in the detector or interactions in the ex- 
ternal surrounding material. The estimated size of the 
intrinsic v e and gamma backgrounds are tied to Mini- 
BooNE event measurements and uncertainties due to 
these constraints are included in the analysis. The in- 
trinsic v e /v e background from muon decay is directly re- 
lated to the large sample of observed v^jv^ events since 
these events constrain the muons that decay in the 50 
m decay region. (The v^jv^ CCQE data sample, in the 
200 < E® E < 1900 MeV energy range, includes 115,467 
and 50,456 neutrino and antineutrino events, respectively 
.) This constraint is accomplished using a joint fit of the 
observed jv^ events and the v e /v e events assuming 
that there are no substantial v^ jv^ disappearance os- 
cillations. The other intrinsic background v e component 
from K-decay is constrained by fits to kaon production 
data and the recent SciBooNE measurements [H] . Other 
backgrounds from mis-identified or z? M [20l [21] events 
are also constrained by the observed CCQE sample. The 
gamma background from NC tt° production mainly from 
A decay or A — > V7 radiative decay [22] is constrained 
by the associated large two-gamma data sample (mainly 
from A production) observed in the MiniBooNE data 
|23j . In effect, an in-situ NC 7r° rate is measured and 
applied to the analysis. Single-gamma backgrounds from 
external neutrino interactions ( "dirt" backgrounds) are 
estimated using topological and spatial cuts to isolate 
these events whose vertex is near the edge of the detec- 
tor and point towards the detector center [3] . 

Systematic uncertainties are determined by consider- 
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FIG. 1: The antineutrino mode (top) and neutrino mode (bot- 
tom) E® E distributions for v e CCQE data (points with sta- 
tistical errors) and background (histogram with systematic 
errors) . 

ing the predicted effects on the u^, v^, v ei and v e CCQE 
rate from variations of parameters. These include uncer- 
tainties in the neutrino and antineutrino flux estimates, 
uncertainties in neutrino cross sections, most of which 
are determined by in-situ cross-section measurements at 
MiniBooNE [20, 23], uncertainties due to nuclear effects, 
and uncertainties in detector modeling and reconstruc- 
tion. A covariance matrix in bins of E® E is constructed 
by considering the variation from each source of system- 
atic uncertainty on the v e and D e CCQE signal, back- 
ground, and and 9^ CCQE prediction as a function of 
E® E . This matrix includes correlations between any of 
the v e and v e CCQE signal and background and and 
CCQE samples, and is used in the x 2 calculation of 
the oscillation fits. 

Fig. [l] (top) shows the E® E distribution for v e CCQE 
data and background in antineutrino mode over the full 
available energy range. Each bin of reconstructed E® E 
corresponds to a distribution of "true" generated neu- 
trino energies, which can overlap adjacent bins. In an- 
tineutrino mode, a total of 478 data events pass the 
P e event selection requirements with 200 < E® E < 
1250 MeV, compared to a background expectation of 
399.6 ± 20.0(sfaf.) ± 20.3(syst.) events. For assessing the 
probability that the expectation fluctuates up to this 478 
observed value, the excess is then 78.4 ± 28.5 events or 
a 2.8(7 effect. Fig. [2] (top) shows the event excess as a 
function of E® E in antineutrino mode. 

Many checks have been performed on the data, includ- 
ing beam and detector stability checks that show that 
the neutrino event rates are stable to < 2% and that 
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FIG. 2: The antineutrino mode (top) and neutrino mode (bot- 
tom) event excesses as a function of E$ E . (E rror bars include 
both the statistical and systematic uncertainties.) Also shown 
are the expectations from the best two-neutrino fit for each 
mode and for two example sets of oscillation parameters. 



the detector energy response is stable to < 1% over the 
entire run. In addition, the fractions of neutrino and an- 
tineutrino events are stable over energy and time, and 
the inferred external event rate corrections are similar in 
both neutrino and antineutrino modes. 

The MiniBooNE antineutrino data can be fit to 
a two-neutrino oscillation model, where the probabil- 



ity, P, of i/„ 



oscillations is given by P 



sm z 28 shr(1.27 Am 2 L/E v ), sin" 20 = 4| i7 e4 1 2 1 ?7 M 4 1 2 , and 



Am 2 



The oscillation parame- 



ters are extracted from a combined fit of the observed 
E® E event distributions for muon-like and electron-like 
events. The fit assumes the same oscillation probabil- 
ity for both the right-sign v e and wrong-sign v e , and 



no significant 



or v e disappearance. Using a 



likelihood-ratio technique [4], the confidence level values 
for the fitting statistic, A\ 2 = x 2 {point) — x 2 (best), as 
a function of oscillation parameters, Am 2 and sin 2 28, 
is determined from frequentist, fake data studies. The 
critical values over the oscillation parameter space are 
typically 2.0, the number of fit parameters, but can be 
as a low as 1.0 at small sin 2 29 or large Am 2 . With 
this technique, the best antineutrino oscillation fit for 
200 < E$ E < 3000 MeV occurs at (Am 2 , sin 2 26) = 
(0.043 eV 2 , 0.88) but there is little change in probabil- 
ity in a broad region up to (Am 2 , sin 2 28) = (0.8 eV 2 , 
0.004) as shown in Fig. [3] (top). In the neutrino oscilla- 
tion energy range of 200 < E$ E < 1250 MeV, the x 2 /ndf 
for the above antineutrino-mode best-fit point is 5.0/7.0 
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with a probability of 66%. The background-only fit has 
a x 2 -probability of 0.5% relative to the best oscillation 
fit and a x 2 / n df = 16.6/8.9 with a probability of 5.4%. 
Fig. [3] (top) shows the MiniBooNE closed confidence level 
(CL) contours for v e and v e appearance oscillations in an- 
tincutrino mode in the 200 < E® E < 3000 MeV energy 
range. The data indicate an oscillation signal region at 
the greater than 99% CL with respect to a no oscillation 
hypothesis, which is consistent with some parts of the 
LSND 99% CL allowed region and consistent with the 
limits from the KARMEN experiment [2"4"j . 

Multinucleon processes and v e and disappearance 
can affect the results of the MiniBooNE oscillation anal- 
ysis. Specifically, nuclear effects associated with neu- 
trino interactions on carbon can affect the reconstruc- 
tion of the neutrino energy, E® E , and the determination 
of the neutrino oscillation parameters [2"5T|2 7j . These ef- 
fects can change the visible energy in the detector and 
the relative energy distribution for the signal and gamma 
backgrounds. These effects are partially removed in this 
analysis since the gamma background is determined from 
direct measurements of NC ir° and dirt backgrounds. 

In order to estimate the possible effects of a 
multinucleon-type model, an oscillation fit was performed 
using event predictions based on the Martini et al. [55] 
model. The prediction was implemented by smearing the 
input neutrino energies as a function of reconstructed 
energy to mimic the behavior of the model. For an esti- 
mate of the effects of disappearance oscillations, a (3+1) 
type model was used. Fits were performed where the ap- 
pearance Am 2 and sin 2 29 app parameters were varied as 
usual but disappearance oscillations were also included 

with \U e4 \ 2 = p^l 2 = \U\ 2 = yan 2 29 app /A and with 
the same Am 2 . This is a disappearance model where all 
four types of neutrinos {y e l^e/ v ^/^^) disappear with the 
same effective sin 2 29di sapP — 4(1 — U 2 )U 2 . A comparison 
of the results for these models versus the nominal Mini- 
BooNE analysis is given in Table|n] Results are presented 
for the best fit with the given prediction model and for a 
test point with Am 2 = 0.5 eV 2 and sin 2 29 = 0.01. The 
difference in x 2 values for the different prediction models 
is < 0.5 units, suggesting that multinucleon or disappear- 
ance effects do not significantly change the oscillation fit 
and null exclusion probabilities. 

Even though the MiniBooNE antineutrino data is a 
direct test of the LSND oscillation hypothesis, the Mini- 
BooNE neutrino-mode data can add additional informa- 
tion, especially for comparisons to various sterile neutrino 
models. The previous MiniBooNE oscillation analysis 
[2] found no evidence for neutrino oscillations in neu- 
trino mode by fitting over the neutrino energy range 
475 < E$ E < 3000 MeV, excluding the low-energy 
range, 200 < E® E < 475 MeV. The reason for excluding 
the low-energy region in the original analysis was based 
on uncertainties for the large gamma background in that 



TABLE II: x 2 values from oscillation fits to the antineutrino- 
mode data for different prediction models. The best fit 
(Am 2 , sin 2 29) values are (0.043 eV 2 , 0.88), (0.059 eV 2 , 
0.64), and (0.177 eV 2 , 0.070) for the nominal, Martini, and 
disappearance models, respectively. The test point \ 2 values 
in the third column are for Am 2 = 0.5 eV 2 and sin 2 28 = 0.01. 
The effective dof values are approximately 6.9 for best fits and 
8.9 for the test points. 





x 2 


values 


Prediction Model 


Best Fit 


Test Pt. 


Nominal V— mode Result 


5.0 


6.2 


Martini et al. [25] Model 


5.5 


6.5 


Model With Disapp. (see text) 


5.4 


6.7 



region. The subsequent work on constraining the low 
energy background and making a more accurate assess- 
ment of the uncertainties now allow the data below 475 
MeV to be used [3]. The neutrino-mode event and ex- 
cess distributions for 6.46 x 10 20 POT are shown in the 
bottom plots of Figs, [I] and [2j respectively. In neu- 
trino mode, a total of 952 events are in the region with 
200 < E® E < 1250 MeV, compared to a background ex- 
pectation of 790.1 ± 28. l(stat.) ± 38.7{syst.) events. This 
corresponds to a neutrino-mode excess of 162.0 ± 47.8 
events with respect to expectation or a 3.4cr excess. 

Two-neutrino oscillation model fits to the MiniBooNE 
neutrino-mode data do show indications of oscilla- 
tions as shown in Fig. [3] (bottom). In contrast to 
the antineutrino-mode results, the MiniBooNE favored 
neutrino-mode region has only small overlap with the 
LSND region and may indicate that the compatibility 
between the two is low in a simple two-neutrino model. 
The best neutrino oscillation fit occurs at (Am 2 , sin 2 20) 
— (3.14 eV 2 , 0.002). In the neutrino oscillation en- 
ergy range of 200 < E$ E < 1250 MeV, the x 2 /ndf for 
the best-fit point is 13.2/6.8 with a fairly small prob- 
ability of 6.1%, and the background-only fit has a x 2 - 
probability of 2% relative to the best oscillation fit and a 
X 2 /ndf = 22.8/8.8 with a probability of 0.5%. As shown 
111 Fig. [2] (bottom), the poor x /ndf for the neutrino- 
mode best fit is due to the data being higher than the 
expectation at low energy and lower at high energy. This 
may be due to the limitation of the simple two-neutrino 
model if the excess is due to oscillations or to some 
anomalous background at low energy if the excess is re- 
lated to backgrounds. 

In summary, the MiniBooNE experiment observes a to- 
tal event excess in antineutrino mode running of 78.4 ± 
28.5 events (2.8cr) in the energy range 200 < E® E < 
1250 MeV. The allowed regions from a two-neutrino fit 
to the data, shown in Fig. [3] (top), are consistent with 
v ^ — > v e oscillations in the 0.01 to 1 eV 2 Am 2 range and 
have some overlap with the allowed region reported by 
the LSND experiment pQ. All of the major backgrounds 
are constrained by in-situ event measurements so non- 
oscillation explanations would need to invoke new anoma- 
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FIG. 3: MiniBooNE allowed regions in antineutrino mode 
(top) and neutrino mode (bottom) for events with E® E > 200 
MeV within a two-neutrino oscillation model. Also shown are 
the ICARUS [28j and KARMEN [24] appearance limits for 
neutrinos and antineutrinos, respectively. The shaded areas 
show the 90% and 99% CL. LSND z? M -> i? e allowed regions. 
The black stars show the MiniBooNE best fit points, while 
the circles show the example values used in Fig. [2] 

lous background processes. The neutrino mode running 
also shows an excess of 162.0±47.8 events (3.4<r), but the 
energy distribution of the excess is marginally compatible 
with a simple two neutrino oscillation formalism. While 
this incompatibility might be explained by unexpected 
systematic uncertainties and backgrounds, expanded os- 
cillation models with several sterile neutrinos can reduce 
the discrepancy by allowing for CP violating effects. On 
the other hand, global fits [T^] with these expanded mod- 
els show some incompatibility with the current upper lim- 
its on electron and muon neutrino disappearance that will 
need new data and studies to resolve. 
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